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1 .  ABSTRACT 


The  linear  and  non-linear  random  modal  interactions  of  a  two  degree-of- 
freedom  aeroelastic  structure  are  examined  by  using  the  Fokker-Planck 
equation  approach.  A  general  differential  equation  describing  the  evolu¬ 
tion  of  the  response  moments  is  derived  for  any  moment  order.  For  the  case 
of  linear  modal  interaction  this  differential  equation  is  found  to  consti¬ 
tute  a  closed  set  of  moment  equations.  The  stationary  response  is  deter¬ 
mined  for  various  system  parameters.  It  is  found  that  the  linear 
interaction  results  in  a  suppression  of  one  mode  when  the  uncoupled  fre¬ 
quencies  of  the  structure  are  close  to  each  other.  For  the  case  of  non¬ 
linear  modal  (known  as  autoparametric)  interaction  the  differential 
equation  of  the  response  moments  forms  an  infinite  coupled  set  of  equations 

which  are  closed  via  two  closure  schemes.  These  are  the  Gaussian  and 
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non-Gaussian  closure  schemes/  The  Gaussian  closure  scheme  requires  14 
coupled  differential  equations  in  the  first  and  second  order  moments,  while 
the  non-Gaussian  closure  leads  to  69  differential  equations  in  the  first 
four  orders  of  response  moments.  The  two  sets  are  solved  by  numerical 
integration.  Both  solutions  exhibit  an  energy  exchange  between  the  two 
modes  in  the  neighborhood  of  the  internal  resonance  condition.  The 
Gauss. inn  closure  solution  gives  a  quasi-stationary  response  in  the  form  of 
fluctuations  between  two  limits.  However,  the  non-Gaussian  closure  solu¬ 
tion  results  in  a  strict  stationary  response.  The  influence  of  random  time 
fluctuations  in  the  system  damping  and  stiffness  coefficients  is  also  exa¬ 
mined.  It  is  found  that  the  damping  variation  has  very  small  effect  on  the 
response  characteristics,  while  the  stiffness  variation  shows  a  pronounced 
effect  on  the  response  mean  squares  for  both  solutions. 
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2 .  RESEARCH  OBJECTIVES 

2.1  BACKGROUND 

The  dynamic  behavior  of  aeroelastic  structures  is  of  main  concern  to 
aeronautical  engineers  who  are  involved  in  the  design  and  reliability  of 
aerospace  structures.  These  structures  are  usually  subjected  to  aerodyna¬ 
mic  forces  which  interact  with  inertia  and  elastic  forces.  The  interaction 
of  these  forces  may  give  rise  to  a  number  of  aeroelastic  phenomena.  For 
example,  the  classical  flutter  (known  also  as  self-excited  oscillation)  can 
occur  due  to  a  linear  interaction  of  these  three  forces.  Classical  flutter 
may  also  involve  the  coupling  of  two  or  more  degrees  of  freedom.  However, 
the  linear  mathematical  modeling  fails  to  predict  a  number  of  observed 
dynamic  characteristics  such  as  amplitude  jump,  limit  cycles,  parametric 
instability,  internal  resonance,  and  saturation  phenomenon.  These  complex 
dynamic  characteristics  owe  their  origin  to  the  inherent  non-linearity  of 
the  structure. 

The  amplitude  jump,  limit  cycles,  and  parametric  instability  a”e  common 
features  of  non-linear  single-  and  multi-degree-of-f reedom  systems. 
Parametric  instability  takes  place  when  the  external  excitation  appears  as 
a  coefficient  in  the  homogeneous  part  of  the  differential  equation  of 
motion.  It  occurs  when  the  excitation  frequency  is  twice  (or  multiple)  of 
the  system  natural  frequency.  Internal  resonance  and  saturation  phenomenon 
may  occur  only  in  non-linear  dynamic  systems  with  more  than  one  degree-of- 
freedom.  Internal  resonance  implies  the  existence  of  a  linear  relationship 
between  the  structure  natural  frequencies  and  results  in  a  non-linear 
interaction  of  the  normal  modes  in  a  form  of  energy  exchange.  Under  exter¬ 
nal  excitation,  the  mode  which  is  directly  excited  exhibits  in  the 
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beginning,  the  same  features  of  a  single  degree-of-f reeaora  system  response 
and  all  other  modes  remain  dormant.  As  the  excitation  amplitude  reaches  a 
certain  critical  level,  the  other  modes  become  unstable  and  the  originally 
excited  mode  reaches  an  upper  bound.  In  this  case  this  mode  is  said  to  be 
saturated  and  the  energy  "spills  over”  into  other  modes.  The  non-linear 
modal  interaction  is  referred  in  the  literature  as  autopararaetric  interac¬ 
tion  since  one  mode  acts  as  a  parametric  excitation  to  another  mode.  The 
non-linearities  in  any  structural  component  may  arise  from  two  main  sources 

1.  Geometric  non-linearities  due  to  large  deformations  such  as  large 
curvature,  end  shortening,  and  inertias  due  to  the  presence  of  con¬ 
centrated  or  distributed  masses. 

2.  Material  properties  which  exhibit  non-linear  or  multi-valued 
stress-strain  relationships.  The  values  of  the  material  properties 
may  also  experience  a  certain  degree  of  uncertainty  due  to  material 
heterogeneity  or  random  temperature  fluctuations.  Thus  the  stiff¬ 
ness  and  damping  coefficients  become  random  variables  or  random 
processes . 

It  is  clear  that  the  aeroe las t ic ian  must  include  the  inherent  structural 
non-linearities  in  his  predictive  raodeis  in  order  to  understand  the  origin 
of  any  unusual  structure  behavior  under  various  types  of  aerodynamic 
loading.  Under  deterministic  unsteady  aerodynamic  forces  these  phenomena 
can  be  predicted  by  one  of  the  standard  techniques  of  non-linear  differen¬ 
tial  equations.  However,  aerospace  structures  are  usually  subjected  to 
turbulent  air  flow,  and  the  analyst  is  encountered  with  aerodynamic  loads 
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which  are  random  in  nature.  These  loads  vary  in  a  highly  irregular  fashion 
and  can  be  described  in  terms  of  statistical  quantities  such  as  means,  mean 
square,  autocorrelation  functions,  and  spectral  density  functions.  The 
dynamic  analysis  of  non-linear  aeroelastic  structures  under  random  loading 
is  not  a  simple  task,  and  it  requires  an  advanced  background  of  probabi¬ 
listic  theory  and  stochastic  differential  equations. 

2.2  MAIN  OBJECTIVES 

In  an  effort  to  understand  the  dynamic  behavior  of  non-linear 
aeroelastic  structures  under  random  excitations  a  research  program  con¬ 
sisting  of  analytical  and  experimental  investigations  is  currently  sup¬ 
ported  by  a  grant  from  the  Air  Force  Office  of  Scientific  Research  (AFOSR) . 
Two  and  three  degree-of-f reedom  systems  possessing  internal  resonance 
are  considered.  The  experimental  investigation  will  demonstrate  the 
existence  of  non-linear  phenomena  and  will  provide  guidelines  for  the  vali¬ 
dity  of  the  theoretical  analysis.  Based  on  the  original  proposal  (February 
1983)  and  its  amendment  (July  1984)  three  main  phases  were  outlined.  These 
are : 

Phase  I:  Investigation  of  the  effects  of  structural  non-linearities 

in  the  neighborhood  of  internal  resonance  conditions  when 
the  structure  is  subjected  to  random  aerodynamic  loading. 

Phase  II:  Investigation  of  the  effects  of  damping  and  stiffness 

uncertainties  in  the  absence  of  internal  resonance  con¬ 
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Phase  III:  Investigation  of  the  combined  effects  of  damping  and  stiff¬ 
ness  uncertainties  in  the  presence  of  internal  resonance 
conditions . 

The  experimental  investigation  will  be  carried  out  on  two  models  emu¬ 
lating  the  analytical  models.  The  models  will  be  excited  by  a  medium  size 
electrodynamic  shaker  of  1200  lb  maximum  thrust  through  a  GenRand  random 
noise  generator.  The  excitation  will  be  Gaussian  filtered  wide  band  pro¬ 
cess  whose  band  covers  a  frequency  range  greater  than  the  normal  mode  fre¬ 
quencies  considered  in  the  analytical  models.  The  excitation  frequency 
band  will  be  adjusted  such  that  higher  modes  will  not  be  excited. 

The  excitation  and  response  processes  will  be  measured  and  recorded 
simultaneously  on  a  magnetic  tape  recorder.  The  mean  squares  and  probabi¬ 
lity  density  of  the  response  will  be  estimated  for  various  values  of  inter¬ 
nal  resonance  detuning.  These  results  will  be  very  valuable  in 
demonstrating  how  the  normal  modes  are  interacting  under  random  excitation. 
In  addition,  the  measured  probability  density  will  be  inspected  for  its  non¬ 
normality  when  the  excitation  is  Gaussian. 

2.3  SUMMARY  OF  MAIN  RESULTS 

The  linear  and  non-linear  random  modal  interactions  of  a  two  degree-of- 
freedom  aeroelastic  structure  are  examined  by  using  the  Fokker  -Planck 
equation  approach.  A  general  differential  equation  describing  the  evolution 
of  the  response  moments  is  derived  for  any  moment  order.  For  the  case  of 
linear  modal  interaction  this  differential  equation  is  found  to  constitute 
a  closed  set  of  moment  equations.  The  stationary  response  is  determined 
for  various  system  parameters.  It  is  found  that  the  linear  interaction 
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results  in  a  suppression  of  one  mode  when  the  uncoupled  frequencies  of  the 
structure  are  close  to  each  other.  For  the  case  of  non-linear  modal  (known 
as  autoparametric)  interaction  the  differential  equation  of  the  response 
moments  forms  an  infinite  coupled  set  of  equations  which  are  closed  via  two 
closure  schemes.  These  are  the  Gaussian  and  non-Gaussian  closure  schemes. 
The  Gaussian  closure  scheme  requires  14  coupled  differential  equations  in 
the  first  and  second  order  moments,  while  the  non-Gaussian  closure  leads  to 
69  differential  equations  in  the  first  four  orders  of  response  moments. 

The  two  sets  are  solved  by  numerical  integration.  Both  solutions  exhibit 
an  energy  exchange  between  the  two  modes  in  the  neighborhood  of  the  inter¬ 
nal  resonance  condition.  The  Gaussian  closure  solution  gives  a  quasi¬ 
stationary  response  in  the  form  of  fluctuations  between  two  limits. 

However,  the  non-Gaussian  closure  solution  results  in  a  strict  stationary 
response.  The  influence  of  random  time  fluctuations  in  the  system  damping 
and  stiffness  coefficients  is  also  examined.  It  is  found  that  the  damping 
variation  has  very  small  effect  on  the  response  characteristics,  while  the 
stiffness  variation  shows  a  pronounced  effect  on  the  response  mean  squares 


for  both  solutions. 
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CHAPTER  I 


BASIC  MODEL  AND  EQUATIONS  OF  MOTION 


1.1  Two-Degree-of-Freedom  System 

Figure  1.1  shows  a  schematic  diagram  of  a  two-degree-of-f reedom 
aeroelastic  structural  model  which  represents  an  aircraft  wing  with  exter¬ 
nal  store.  The  model  consists  of  two  coupled  beams  of  stiffnesses,  k^  and 
and  two  tip  masses,  m^  and  .  When  the  horizontal  beam  is  subjected 
to  a  random  motion  £(t),  the  two  beams  will  move  as  shown  in  Fig.  1.1 
with  tip  deflections  q^  and  ^ ,  respectively.  By  applying  Lagrange's 

•k 

equation,  and  including  the  effects  of  the  axial  shortening  motions,  q^ , 

* 

and  q^  in  the  kinetic  energy  expression,  the  equations  of  motion  are  given 
in  terms  of  the  generalized  coordinates  q^  and  q2 : 
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Non-lineartties  up  to  cubic  order  are  retained  in  ^  ^  and  4'2‘  It  seen 

that  the  equations  of  motion  contain  linear  dynamic  coupling  since 

mij  *  I  *  J*  In  addition,  the  support  motion  acceleration  £(t)  enters 

the  equations  of  motion  as  a  non-homogeneous  forced  excitation  l  ( t )  and 

as  a  parametric  excitation  as  given  by  the  terms  £  ^  (t)q^.  The  functions 

and  \l>2  include  all  nonlinear  inertia  terms.  These  can  be  subdivided 

into  two  groups.  The  first  group  represents  non-linear  inertia  of  the  same 
2" 

mode  such  as  q^q^,  in  the  first  equation.  The  second  group  represents  auto- 
parametric  coupling  (or  non-linear  interaction)  such  as  q^’c^  in  the  first 
equation.  Here  the  acceleration  (which  is  an  implicit  function  of  time) 
acts  as  a  parametric  excitation  to  q^  motion. 

In  deriving  equations  (1.1)  it  was  assumed  that  both  damping  coef¬ 
ficients  C.^  and  beam  stiffness  are  constant  coefficients.  In  Chapter 

4  these  coefficients  will  be  subject  to  random  fluctuations. 

In  Chapter  2,  the  linear  modal  Interaction  under  random  excitation  will 
be  examined  by  dropping  ^  and  In  Chapter  3,  the  non-linear  modal 

interaction  will  be  investigated  when  the  structure  is  tuned  to  the  inter¬ 
nal  resonance  condition  =  2w^  (where  and  u>2  are  two  eigenvalues  of 
the  structure).  Chapter  4  will  treat  the  influence  of  random  fluctuations 
of  the  system  parameters. 


CHAPTER  II 


LINEAR  ANALYSIS 


11. 1  Introduction 

The  normal  mode  interaction  of  the  two-degree-of-f reedom  system  shown 
in  Fig.  1.1  will  be  examined  within  the  framework  of  the  linear  theory  of 
random  vibration.  In  other  words,  the  system  response  represented  by 
equations  (1.1)  will  be  analyzed  after  setting  the  non-linear  functions 
^  and  ip 2  to  zero.  The  normal  mode  frequencies  and  mode  shapes  will  be 
obtained  first,  and  then  the  equations  of  motion  will  be  transformed  into 
normal  coordinates.  These  in  turn  will  be  written  in  the  Ito  type 
stochastic  differential  equations  [1,2]. 

The  differential  equations  of  the  response  statistical  moments  will  be 
derived  by  using  the  Fokker-Planck  equation  [1].  The  Fokker-Planck 
equation  is  a  partial  differential  equation  which  describes  the  evolution  of 
the  system  response  probability  density  with  respect  to  time  and  response 
coordinates.  In  view  of  the  linearity  of  the  system  equations  of  motion, 
the  moment  equations  will  form  a  closed  set  which  can  be  solved  numerically 
by  one  -.if  the  standard  techniques.  The  mean  square  of  the  response  moment 
will  be  examined  for  a  wide  range  of  system  parameters.  The  influence  of 
random  parametric  coefficients  will  also  be  examined. 

11. 2  Normal  Mode  Analysis 

By  dropping  the  non-linear  function  <|^  and  ^  from  the  equations  (1.1), 


the  system  becomes  linear  with  random  coefficients: 
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of  the  mass  parameter  p . 


Equations  (II. 1)  can  be  transformed  into  normal  coordinates  via  the 


The  following  non-dimensional  parameters  are  introduced 


Y  =  !i 

1  o 
ql 


v  -I* 
2 


(II. 6) 


where  is  the  root  mean  square  response  of  the  system  when  it  is  reduced 

into  a  single  degree-of-f reedom  represented  by  the  horizontal  beam  with  end 

2 

mass  (m^  +  m2).  Thus,  E[q°]  is  the  mean  square  response  of  the  single- 
degree-of-f reedom  system. 


..  O  |  rt  O,  2  O  ,  ,y  , 

ql +  i’lla)llql  +alllql =  ” w(c) 


(II. 7) 


The  random  excitation  acceleration  £(t)  has  been  replaced  by  the  zero 
mean  Gaussian  white  noise  process  W(t)  whose  autocorrelation  function  is 
defined  by  the  relation 


R[A  t  ]  =  E[W(t)W(t  +  At)]  =  2D  6  (At) 


(II. 8) 


where  2D  is  the  spectral  density  of  the  random  process  W(t)  and  6 (  )  is  the 
Dirac  delta  function. 

The  mean-square  response  of  system  (II. 7)  is  given  by  the  well-known 
solution  [3] : 


E[q°j  = 


2^U  «u 


Thus,  we  may  select  q^  to  be  the  root  mean  square  of  e[ q°] ;  i.e., 


(II. 9) 


2c  11  “ll 


E[  ]  denotes  expectation 


The  equations  of  motion  in  terms  of  the  non-dimensional  normal  coor¬ 


dinates  Y.  are 
1 

Y1  +  2ciYi  +  Y1  =  riEl  +  £(E2Y1  +  E3Y2)]m(t) 

2  (II. 10) 

Y2  +  2rC2Y2  +  r  y2  =  tGi  +  c(G2Y2  +  G3Y2)]w(t) 


These  are  two  non-homogeneous  stochastic  differential  equations  which  are 
coupled  through  the  parametric  terms  Y^W(x )  and  Y2W(t).  The  prime  denotes 
differentiation  with  respect  to  the  non-dimensional  time  x.  The  coef¬ 
ficients  E^,G^  are  defined  in  Appendix  A. 


II.  3  Response  Moment  Equations 

The  Fokker-Plank  equation  or  the  Ito  stochastic  calculus  can  be  used  to 
generate  the  differential  equations  of  the  response  moment  as  outlined  by 
Ibrahim  [ A ] . 

Introducing  the  state  variable  transformation 

fY0  fx0 


equations  (II. 10)  can  be  written  in  the  state  Markov  form 

xj  -  x3 

*2  =  X4 

(11.12) 

I 

x3  =  -  Xj  -  2^jX3  +  (F. ,  +  cF.2Xl  +eE3X2)  W(t) 

'  2 

X  =  -  r  X  -  2 r^rX,  +  (G,  +  tGnX,  +  eG^XJ  W(t  ) 


The  non-stat ionary  Fokker-Planck  equation  of  the  evolution  of  the  pro¬ 


bability  density  of  the  response  vector  X  is 


(X,T)  =  -  X-  (  X->P(X,t  )}  -  £ -  {x,p(x,t)} 


ax  1  3l 


ax,  1  4 


{  [-X1  -  2?1X3]p(X,t)} 


~  {[-  ^2X2  -  2c2rX4]P(X,T)} 

4 


+  D  ~  {[e2  +  2eE1E2X1  +  2eE1E3X2  +  e  E2X2 


+  2eZE2E3X1X2  +  e2E3X2]  P(X,t  )} 


(11.13) 


+  20  alO-xT  UVi  +  £(E1G2  +  E2G1)X1  +  £(E1G3  +  E3G1)X2 
3  4 


+  e2E2G2X^  +  e2(E2G3  +  E3G2)X1X2  +  e  ^E^X*]  p(X ,x  )} 


2j-.  ~  v2i 


+  D  i-g-  {[G2  +  2eGLG2X1  +  2eG1G3X2  +  e  2G2X2 

8  X  # 

4 


4  2£*G2G3XlX2  +  £2G2x2]p(X,t)} 


The  following  notation  for  the  response  joint  moment 


00 

■iiw  -////  Sllx2XX  Pff.T)dX1dX2dX3dX4 


(11.14) 


will  be  adopted. 
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Premultiplying  both  sides  of  equation  (11.13)  by  x|x^X^X^  and 
integrating  both  sides  over  the  whole  space  -»<X<"  gives  the  following 
general  moment  differential  equation: 


mijk£  imi-l,  j,k+l,£  +  jmi,  j-l.k.A+l 


k(mi+l,j,k-l,A  +  2^lmi,  j,k,£) 


-A(r  “i.J+l.M-l  +  *2r"i,j,M> 


+  k(k-l)D(E^.jj>k_2>jl  +  2eE1E2«i+1>jtk.2 


+  2eE1E3m.^+1^k_2^  +  ^E2mi+2 ,  j  ,k-2 ,1 


'  E2E3mi+l,  j+l,k-2,£  +e  E3mi ,  j+2  ,k-2  ,lj 


+  2k£D[ElGlmi,j,k-l,£-l  +  <E1G2  +  E2Gl)mi+l, 

+  £(E1G3  +  E3Gl)mi,  j+l.k-l.t-1  +  e2E2G2mi+2,j,k-l,£-l 
+  e  (E2G3  +  E3G2)m.+1(  j+1>k_3  +  e2E3G3mi,  j+2,k-l,£-l 


+  e 


E3G3rai, j+2,k-l,£-l] 


+  U*'1)D  j,--2  +  2£GlG2Vl,j,k^2 


+  2eG1G3rai>  j+l>kj£_2  +  e2G2mi+2,  j,k,*-2 


+  2e  G2G3mi+i  t  j+\  ,k,£-2  +e  G3mi ,  j+2  ,k,£ -2^ 


(11.15) 
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An  inspection  of  equation  (11.15)  shows  that  the  moment  equations  of  any 
order  are  consistent,  that  is,  the  moment  equations  of  order  n  are  not 
coupled  with  moments  of  order  greater  than  n  and  they  can  be  solved  analy¬ 
tically. 

The  first  and  second  order  response  moments  are  of  main  concern  in  ran¬ 
dom  response  analysis.  The  first  and  second  order  moment  equations  are 
generated  from  (11.15)  and  the  following  stationary  solution  is  obtained: 

m1000  =  m0100  =  m0010  =  m0001  =  0 


"1010  “  “0101  ~  0 


m0002  r  ni0200 


(11.16) 


m0110  m1001 


1  -  r 

m1001  =  2(c L  +  r?  )  “llOO 


12000 ’  m0200 ’  m1100  are  8iven  ^  solution  of  these  three  equations: 


d!  d2  d3 


el  e2  e3 


fl  f2  f3 


1100  )  (  f4 


(11.17) 


where  d^,  e^,  f ^  are  defined  in  Appendix  B. 

The  solution  of  equations  (11.17)  for  ji^qqq  and  m^gg  is  sh°wn  in 
Fig.  II. 3  as  a  function  of  the  frequency  rat  io  u) ^  for  various  values 
of  the  mass  parameter  p.  The  response  mean  square  (in  terms  of  the  genera¬ 
lized  coordinates)  is  obtained  through  the  inverse  modal  transformation  and 


is  plotted  in  Fig.  II. 4.  Figures  II. 3  and  II. 4  reveal  that  the  mean  square 


response  exhibits  a  strong  interaction  in  the  neighborhood  of  the  frequency 
ratio  1 1  =  * *  The  degree  °f  interaction  depends  on  the  mass  parameter 

p  which  measures  the  degree  of  linear  dynamic  coupling.  The  interaction  of 
the  two  modes  exhibits  the  linear  vibration  absorber  characteristics  known 
in  the  deterministic  vibration.  It  is  also  seen  that  as  the  mass  parameter 
p  diminishes,  the  absorbing  effect  increases.  For  modest  values  of  p  <1, 
the  first  mode  is  suppressed  and  the  second  mode  Y2  reaches  its  peak 
value.  Furthermore,  the  response  of  the  main  horizontal  beam  is  that  of  a 
single-degree-of-freedom  system  outside  the  frequency  ratio  u22^wll  =  1  - 
0(e),  and  the  range  of  this  frequency  ratio  is  influenced  by  the  mass  ratio 
parameter  p  as  shown  in  Fig.  II. 4. 

In  order  to  examine  the  influence  of  parametric  random  excitation  on 
the  system  response,  the  stationary  solution  was  obtained  by  setting 

i  =  1,2  to  zero  from  the  moment  equations.  Figures  II. 5  and  II. 6 
show  the  mean  square  response  in  terms  of  the  normal  coordinates  and 
generalized  coordinates,  respectively.  It  is  seen  that  the  parametric 
excitation  has  very  negligible  effect  on  the  overall  response  level. 
However,  it  is  the  system  stochastic  stability  which  is  governed  by  the 
parametric  excitation  [4] . 


external  excitation 


CHAPTER  III 


NON-LINEAR  ANALYSIS 


III.l  Introduction 

The  linear  analysis  of  the  coupled  aeroelastic  structure  examined  in 
Chapter  2  may  be  sufficient  to  describe  the  dynamic  behavior  of  the  system 
as  long  as  its  motion  is  very  close  to  the  static  equilibrium  position. 
However,  the  structure  may  experience  new  types  of  modal  interactions  when 
the  non-linear  effects  become  significant.  Of  particular  importance  is  the 
non-linear  coupling  of  the  system  normal  modes.  It  is  well  known  in  the 
deterministic  theory  of  non-linear  oscillations  [5]  that  this  type  of  non¬ 
linearity  may  give  rise  to  the  dynamic  instability  known  as  "internal 
resonance."  Methods  for  studying  non-linear  effects  in  random  vibration 
problems  involve  a  number  of  difficulties  in  determining  the  response  dyna¬ 
mic  characteristics  [4] .  These  difficulties  include  the  solution  of  the 
system  Fokker-Planck  equation  and  the  problem  of  infinite  coupled  moment 
equations  of  the  response  coordinates. 

Over  the  past  few  decades  a  number  of  closure  schemes  have  been  deve¬ 
loped  [6] .  In  the  present  analysis  two  closure  schemes  will  be  proposed  to 
close  the  me  non:  dynamic  equations  of  the  response  coordinates.  The  first 
scheme  is  based  on  the  assumption  that  the  response  process  is  Gaussian 
distributed  under  Gaussian  excitation.  However,  the  application  of  this 
method  to  non-linear  systems  is  not  mathematically  justified  since  the 
response  is  non-Gaussian.  Therefore,  it  is  assumed  that  the  excitation  is 
of  small  intensity  and  the  response  will  not  depart  significantly  from  nor¬ 


mality.  The  response  process  can  then  be  described  adequately  in  terms  of 
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the  first  and  second  order  cumulants  (semi-invariants),  and  all  higher 


order  cumulants  vanish.  Cumulants  are  statistical  functions  which  can  be 
related  to  the  moments.  For  example,  first  and  second  order  cumulants  are 
equivalent  to  the  mean  and  variance  of  the  response,  respectively. 

The  second  method  is  more  accurate  and  includes  the  effects  of  the  non¬ 
normality  of  the  response.  Unlike  Gaussian-closure  methods,  cumulants  of 
order  higher  than  two  will  not  vanish  and  will  give  a  measure  of  the 
deviation  of  the  response  from  normality.  To  the  first  order  approximation 
third  and  fourth  order  cumulants  will  be  significant  to  account  for  the 
response  non-normality  while  the  fifth  order  cumulant  will  contribute  less 
and  can  be  neglected.  Consequently,  one  can  express  the  fifth  order  joint 
moments  in  terms  of  lower  order  moments.  Hitherto,  the  present  approach 
has  not  been  applied  to  multi-degree-of-f reedora  non-linear  systems. 

It  is  expected  that  new  features  of  the  response  characteristics  may  be 
obtained  due  to  the  system's  inherent  non-linearity.  The  influence  of  the 
internal  resonance  condition  will  be  examined  for  various  values  of  the 
system  parameters. 


III. 2  Theoretical  Analysis 


The  non-linear  equations  of  motion  of  the 
1.1,  are 


1-  2 

i  ! 

..  / 

mL  +  m2(l  +  2.25(ji.)  ) 

1-5-2  if 
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system  (2.5),  shown  in  Fig. 


a 


(in.i) 


where 


n  n  £■  "  .  n  /  c  '  2  . 2.  1.2  f  "  ,  .  Z',  ,  U  .  0 

*1  =  °*9  7T  qlql  +  °-45  -I  ql  +  f—  (  q2q2  +  q2^  +  rr  qlq2 


.2i  ,  Ch3 

1 


+  r^Kq2  +  qiq2^ 


0.3  "  1.2.2  1.2 

*2  ‘  q-  V.  -  *T  "i  +  q-  Vi 


k. 

1  7^ 

l 


i  =  1,2 


Equations  (III.l)  can  be  written  In  terms  of  normal  coordinates  y.  via  the 
transformation 


{q}  =  [R]  { y} 


(III. 2) 


where  modal  matrix  [R]  was  defined  in  section  II. 2. 

Premultiplying  equation  (III.l)  by  [R]  1  [m]  * ,  where  [m]  is  the  mass 
matrix  and  using  transformation  (III. 2),  the  equations  of  motion  take  the 
form 


(III. 3) 
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where  prime  denotes  differentiation  with  respect  to  the  time  parameter 
t  =  u)-^t  and  r  =  A  linear  viscous  damping  has  been  introduced  to 

account  for  energy  dissipation  of  the  system.  The  following  non- 
dimensional  parameters  have  been  used: 


(III. 4) 


(y1>y2}  =  { y1*y2^ /ql 

q°  is  the  response  mean  square  of  the  system,  which  was  defined  in  section 
II. 2,  and  the  non-linear  functions  and  ^  are 


h  -  VlY'l  +  VlV2  +  %Vl  *  a7Y2Y2  +  Vi2  +  a9YlY2  +  a10Y^ 


(III. 5) 

*2  ‘  b4 Y 1 Y 1  +  VlY2  +  b6Y2Yl  +  bJY2Y2  +  V?  +  VA  +  be*2 

where  the  coefficients  a^  and  are  given  in  Appendix  C.  Two  types  of 

non-linearity  are  embodied  in  The  first  forms  the  non-linear  terms 

••  »2 

of  the  same  mode  such  as  Y^Y^  or  Y^  in  the  first  equations  of  mode  1,  and 
the  second  constitutes  autopararaetric  terms  such  as  Y^Y^  (non-linear 

u  2 

coupling).  Autoparametric  terms  give  rise  to  r  =  —  =  0.5,  (here 
where  the  ord.r  of  modes  is  reversed). 

The  acceleration  £"(t)  is  assumed  to  be  a  Gaussian  wide  band  random 
process  with  zero  mean  and  a  smooth  spectral  density  2D  up  to  some  fre¬ 
quency  which  is  higher  than  any  characteristic  frequency  of  the  system. 

In  order  to  represent  the  response  coordinates  as  a  Markov  process  the 
acceleration  Y^  associated  with  the  non-linear  terras  must  be  removed  by 
successive  elimination.  Having  eliminated  from  the  non-linear  terras  in 
equations  (III. 3)  the  following  coordinate  transformation  is  introduced: 


(III. 6) 

Equations  (III. 3)  can  be  written  in  the  Markov  vector  form: 

Xi  =  X3 
X2  =  X4 

x3  =  ~  x!  "  2^X3  -  a4xj  +  (a6  +  r2a5)  X^  -  r2a?X2 

-  2c1a4X1X3  -  2?2ra5X1X4  -  2?1a6X2X3  -  2^2ra?X2X4  +  agX2 

+  a9X3X4  +  *10*5  -  W  +  A2X!  +  A3X2  +  A4X5  +  A5X1X2  +  A6X5> U<T > 

*4  '  -  'S  -  *lrXt  -  b4X5  -  ("6  +  X1X2  ‘  r%X5  ‘  * lb4XlX3 

2^2rb5XlX4  "  24 1 b6X2X3  "  2c2rb7X2X4  +  b8X3  +  b9X3X4  +  b10X4 

-  (bl  +  B2X1  +  B3X2  +  B4X2  +  63X3X3  +  B6X2)  W(t) 

ar.'  7) 

where  the  coefficients  A^  and  B^  are  given  in  Appendix  D. 

In  the  equations  (III. 3)  the  random  acceleration  £  "(t )  has  been 
replaced  by  the  white  noise  process  W(x),  where  the  Wong-Zakai  [7]  correc¬ 
tion  terra  is  zero.  The  autocorrelation  function  of  W(t )  is  defined  by  the 
relationship 

R[At1  =  E[W(t)W(t  +  At)]  =  2D  6  (Ax  )  (III. 8) 

where  2D  is  the  spectral  density,  and  6(  )  is  the  Dirac  delta  function. 


r.  - 


.V 


28 


L 


A  general  differential  equation  for  all  possible  moments  can  be 
generated  by  using  the  Ito  stochastic  calculus  or  the  Fokker-Planck 
equation  [4].  Following  the  same  procedure  described  in  Section  II. 3  the 
general  differential  equation  of  the  response  joint  moment  is 


mi,j,k,fl  imi-l ,  j ,k+l  ,£  +  'imi,  j-1  ,k,£+l 


+  k[  rai+l,  j,k-l,«.  2?lmi,j,k,S, 
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-  r  a_m. 
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(HI. 9) 


H 


.7,  ^  »  L 


,  .  *  -  WV *\  V.  V,.  ,  „  .  .  .  _ 

J  >  -  ~-a 


+  2d(A1B4  +  A4Bl  +  A2B2Jm.+2>j>k_lt£_1 
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+  2DBlB3"l,1+l,k,t-2  +  *2B1B4  +  B2)"i+2,j,k,t-2 
+  2D(BlB;  +  B2B3)  mi+l ,  j+l  ,k,Jt-2  +  *  B?  +  2B1B6>  “l ,  j+2  .k.t-l' 

(III. 9) 
cont ' d . 

It  is  seen  that  a  moment  equation  of  order  n  =  i+j+k-W.  contains  moments  of 
order  n  and  n+1 ,  thus  forming  an  "infinite  hierarchy  set."  In  order  to 
close  the  moment  equations  the  following  two  closure  schemes  will  be  used. 


III. 3  Gaussian  Closure  Solution 


From  equation  (III. 9)  it  is  possible  to  generate  four  equations  for  the 
first  order  moments  and  ten  equations  for  the  second  order  moments.  These 
equations  are,  however,  coupled  through  third  order  moment  terras. 

Given  the  assumption  that  the  system  non-linearities  are  too  small  to 
the  extent  that  the  response  can  be  regarded  as  "nearly"  Gaussian,  then  the 
fourteen  equations  can  be  closed  by  using  the  Gaussian  cumulant-negloct 
scheme  [4] .  Under  this  assumption  the  cubic  semi-invariants  vanish  and  the 
third  order  moment  terras  can  be  expressed  in  terms  of  lower  order  moments. 


Mv/J  =  EtXiXjXJ  '  l  E[Xil  EC  XjXk]  +  2E[  x±]  e[  Xj]  F.[  Xk]  =  0  (III. 10) 
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where  the  number  over  the  summation  sign  refers  to  the  number  of  terms 
generated  by  the  indicated  expression  without  allowing  permutation  of 
indices . 

The  resulting  closed  14  coupled  non-linear  differential  equations  are 
integrated  numerically  by  using  IMSL  DVERK  routine  (Runge-Kutta-Verner 
fifth  and  sixth  order  numerical  integration  method). 

Figure  III.l  shows  the  time  history  of  the  system  mean  square  response 
in  normal  coordinates  for  internal  resonance  ratio  r  =  0.5,  damping  ratios 
^ l  =  ^  2  =  mass  ratio  m^/m^  =  0.2,  beam  length  ratio  l ^  =  0.6  and 

non-linear  coupling  parameter  e  =  0.02.  After  a  sufficient  period  of  time 
t  =  1000  the  response  fluctuates  between  two  limits,  indicating  that  the 
system  does  not  achieve  a  stationary  response. 

The  effect  of  damping  ratios  and  5 ^  on  the  response  mean  squares  in 
normal  and  generalized  coordinates  is  shown  in  Fig.  III. 2.  It  is  seen  that 
as  the  damping  ratios  decrease,  the  region  of  autoparametric  interaction 
becomes  wider  and  the  peak  of  the  mean  square  response  of  the  two  modes 
increases.  Also,  the  quasi-stationarity  of  the  Gaussian  closure  solutions 
are  manifested  over  a  wider  range  of  r  = 

Figure  III. 3  shows  the  effect  of  the  non-linear  coupling  parameter  z. 
For  very  small  values  of  s  ,  the  effect  of  non-lineuriti in  the.  equations 
of  motion  (III.l)  is  greatly  reduced.  The  mean  square  responses  do  not 
exhibit  any  non-linear  modal  interaction  and  almost  follow  the  linear  solu¬ 
tions.  However,  it  can  be  seen  that  the  system  responds  very  differently 
for  even  minor  increases  in  e .  As  e  increases,  the  interaction  region 
becomes  wider  and  the  autoparametric  interaction  takes  place  in  a  form  of 
energy  exchange  between  the  two  modes. 


The  effect  of  mass  ratio  is  plotted  in  Fig.  III. 4.  As  the  mass  ratio 
is  increased,  the  region  of  autoparametric  interaction  between  the  displa¬ 
cement  mean  square  responses  in  normal  coordinates  shrinks  to  a  narrow 
range  of  internal  tuning  ratio. 

Figure  III. 5  illustrates  the  effect  of  length  ratio  between  the  two 
beams  on  the  mean  squares  of  the  system  response.  As  the  length  ratio  is 
increased,  the  region  of  autoparametric  interaction  shrinks  to  a  narrow 
range  of  internal  tuning  ratio. 


Normal 


III. 5  Gaussian  closure  solution  for  various  values 
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III. 4  Non-Gaussian  Closure  Solution 

In  Che  theory  of  random  processes  [8]  it  is  known  that  any  linear 
operator  on  a  Gaussian  process  results  in  another  Gaussian  process. 
However,  if  the  system  is  non-linear  or  involves  random  coefficients,  then 
the  response  will  not  be  Gaussian  and  the  corresponding  third  and  higher 
order  semi-invariants  will  not  vanish.  These  higher  order  semi-invariants 
(cumulants)  give  a  measure  to  the  non-normality  of  the  response  distribu¬ 
tion.  If  the  response  process  is  assumed  to  be  slightly  deviated  from 
Gaussian,  the  contribution  of  the  higher  order  cumulants  diminishes  as  the 
order  increases.  Under  this  assumption,  a  better  approximation  can  be 
obtained  by  letting  the  fifth  and  higher  order  cumulants  be  zero,  i.e., 


StVjWnJ  '  EKXjW„)  -  I  ElXl)EtXjXkXtxJ 


10  10 
+  2  I  EfxjEfX^EfX^xJ  -  6  I  E[  xj  e[  X.]  E[  xj  e[  Xj 


(III. 11) 


15  10 

+  2  l  EtxjEfx.XjE^XJ  -  l  Efx^Efx^xJ 

+  24E[Xi]E[X.]E[xk]E[xjE[xJ  =  0 


Moment  differential  equations  of  order  up  to  four  will  be  generated 
from  equation  (III. 9).  Sixty-nine  of  these  equations  are  coupled  ough 
the  fifth  order  moment  terras.  These  69  equations  are  closed  by  using  the 
relation  expressed  in  (III. 11). 

These  equations  are  integrated  numerically  by  using  TMSL  DVERK  routine 
(double  precision).  The  time  history  response  of  the  mean  square  displace¬ 


ments  in  normal  coordinates  are  shown  in  Fig.  II. 6  for  internal  resonance 
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ratio  r  =  0.5,  damping  ratios  t;  ^  ~  0*02,  mass  ratio  nij/m^  =  0.2  and 

length  ratio  SL  ^1^  \  ~  0.6.  During  the  transient  period  the  mean  square  of 
the  first  normal  mode  displacement  fluctuates  and  grows  until  it  reaches  a 
peak  value  at  x  =60  and  then  drops  to  a  lower  level  at  x  =  110.  During 
this  transient  period,  the  mean  square  of  the  second  normal  mode  displace¬ 
ment  fluctuates  and  drops  until  it  reaches  its  minimum  value  at  x  =  60, 
then  grows  not  significantly  to  a  peak  value  at  x  =  110.  This  kind  of 
interaction  shows  an  energy  exchange  between  the  two  modes  during  a  tran¬ 
sient  period,  after  which  each  mode  achieves  a  complete  stationary 
response . 

The  stationarity  of  the  solution  is  confirmed  by  solving  numerically 
the  non-linear  algebraic  equations  resulting  from  the  set  of  original  dif¬ 
ferential  equations.  The  numerical  solution  is  achieved  by  using  the  TMSL 
routine  ZSCNT  (Secant  method  for  simultaneous  non-linear  equations).  The 
algebraic  numerical  solution  is  identical  to  the  stationary  solution 
obtained  by  numerical  integration.  However,  for  Gaussian  closure,  the 
algebraic  solution  does  not  converge  for  all  possible  initial  guessing 
values . 

Figure  III. 7  shows  the  effect  of  damping  ,".atios  on  the  autoparametric 
interaction  r  gion.  It  is  seen  that  as  the  damning  ratios  decrease  the 
region  of  autoparametric  integration  between  the  two  modes  broadens  and  the 
difference  between  the  peak  mean  square  responses  of  the  two  modes 
increases,  but  less  than  in  the  Gaussian  closure  solution.  Unlike  the 
Gaussian  closure  solution,  the  stationarity  of  the  non-Gaussian  closure 
solution  is  manifested  over  a  wider  range  of  r  = 


r 
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The  effect  of  the  non-linear  coupling  parameter  e  is  shown  in  Fig.  III. 8. 
The  small  value  of  e,  which  exhibits  no  non-linear  effect  in  the  Gaussian 
closure  solution,  shows  modal  interaction  in  the  non-Gaussian  closure  solu¬ 
tion.  Also,  as  the  non-linear  coupling  parameter  e  increases,  the  region 
of  modal  interaction  widens  and  the  difference  between  the  two  peak  values 
increases . 

The  effect  of  the  mass  ratio  on  the  mean  square  responses  of  the  system 
is  shown  in  Fig.  III. 9.  As  the  mass  ratio  increases,  the  mean  square  of 
the  second  mode  displacement  increases,  accompanied  by  a  suppression  in  the 
mean  square  of  the  first  mode.  Graphs  are  shown  for  both  normal  and 
generalized  coordinates. 

Figure  III. 10  demonstrates  the  effect  of  length  ratio  on  the  mean 
squares  of  the  system  response.  As  the  length  ratio  increases,  the  mean 
square  of  the  second  mode  displacement  decreases.  This  decrease  is  accom¬ 
panied  by  a  noticeable  increase  in  the  mean  square  of  the  first  mode 
displacement . 

The  effect  of  initial  conditions  on  the  system  mean  square  response  is 
then  examined  to  see  if  the  system  possesses  more  than  one  limit  cycle. 

»  t 

When  the  initial  condition  of  E[  (=  is  set  equal  to  0.00001 

and  all  otner  moments  have  zero  initial  conditions,  the  Gaussian  closure 
solution  gives  the  same  quasi-stationary  response  as  those  solutions 
obtained  with  a  different  set  of  initial  conditions  (compare  Figs.  III.l 
and  III. 11).  In  the  case  of  the  non-Gaussian  closure  solutions,  it  is  also 
found  that  the  initial  conditions  have  no  effect  on  the  final  steady  state 
response,  as  shown  in  Figs.  III. 6  and  III. 12. 


The  C.P.U.  time  for  I.B.M.  3033  system  required  to  generate  the 
Gaussian  solution  up  to  x  =  1500  is  24  seconds  while  for  the  non-Gaussian 
closure  solution,  it  is  487  seconds  for  the  same  response  period. 


Fig.  III. 7  Non-Gaussian  closure  solution  for  various  values  of  damping  ratio. 


closure  solution  for  various  values  of  non-linear  coupling  parameter 


Fig.  III.'}  Non-Gaussian  closure  solution  for  various  values  of  mass  ratio 


Ill .5  Conclusions 


The  random  response  of  an  aeroelastic  structure  with  autoparametric 
coupling  has  been  investigated  in  the  neighborhood  of  the  internal  reso¬ 
nance  condition  r  =  =  0.5.  The  response  has  been  determined  via  two 

closure  schemes:  Gaussian  and  non-Gaussian.  In  the  time  domain,  the  tran 
sient  response  in  both  schemes  exhibits  the  well-known  characteristics  of 
autoparametric  interaction  in  a  form  of  energy  exchange  between  the  dis¬ 
placement  mean  squares  of  the  two  modes.  Furthermore,  the  level  of  these 
mean  squares  may  exceed  the  mean  squares  during  the  steady  state  period. 
The  Gaussian  closure  solution  yields  a  quasi-stationary  response  while  the 
non-Gaussian  closure  solution  gives  a  stationary  response. 

The  stationarity  of  the  non-Gaussian  closure  solution  may  be  confirmed 
from  the  analytical  solution  of  Schmidt  [9]  who  employed  the  stochastic 
averaging  method  to  a  non-linear  two-degree-of-f reedom  system.  However, 
Schmidt  could  not  derive  a  closed  form  expression  for  the  normalized 
constant  of  integration  of  the  response  probability  density.  The  present 
investigation  gives  the  dynamicist  some  guidelines  for  using  the  Gaussian 
or  non-Gaussian  closures  in  more  complicated  systems  where  none  of  the  ana 
lytical  approaches  is  applicable. 


CHAPTER  IV 


NON-LINEAR  RANDOM  RESPONSE  IN  THE 
PRESENCE  OF  PARAMETER  UNCERTAINTY 


IV. 1  Introduction 

The  dynamic  behavior  of  an  aeroelastic  non-linear  two-degree-of-f reedom 
system  has  been  examined  in  Chapter  III.  The  random  excitation  appeared  in 
the  equations  of  motion  as  a  non-homogeneous  term  and  as  a  parametric  coef¬ 
ficient.  It  was  assumed  that  both  damping  and  stiffness  properties  of  the 
structure  were  time  independent  (constants).  As  the  structure  oscillates 
these  properties  may  experience  a  certain  degree  of  fluctuation  as  a  result 
of  the  inherent  temperature  variation.  Measurements  taken  in  laboratory 
experiments  often  show  that  dynamic  properties  such  as  damping  and  stiff¬ 
ness  of  vibratory  systems  are  non-repeatable  parameters;  every  experiment 
gives  different  values  for  damping  and  stiffness  of  the  same  system,  and 
the  differences  are  random. 

The  dynamic  response  of  structures  with  randomly  varying  parameters  is 
of  great  practical  interest  to  the  analyst  or  designer  using  lumped  para¬ 
meter  models,  because  the  random  variation  of  some  system  parameters  may 
cause  the  system  respone  to  exceed  design  specifications.  These  random 
uncertainties  can  be  classified  into  two  main  categories  MG'!:  statistical 
and  non-statistical.  To  cite  just  one  example,  statistical  uncertaint res 
can  be  due  to  the  stiffness  or  damping  fluctuations  caused  by  random 
variation  in  material  properties  or  variation  caused  by  manufacturing  and 
assembly  techniques.  Non-statistical  uncertainties  can  be  due  to  the 
approximation  involved  in  the  mathematical  modeling  of  complex  structural 
systems . 
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It  is  very  important  to  distinguish  between  two  different  problems 
encountered  in  mechanical  vibrations  and  aeroelastic  flutter.  These  are 
the  random  response  of  dynamic  systems  to  random  parametric  excitations 
which  appear  as  coefficients  in  the  equations  of  motion,  and  the  random 
response  to  random  external  excitations  when  the  system  parameters  are  not 
precisely  defined  and  represented  in  a  probabilistic  sense.  In  the  former 
case  the  system  equations  of  motion  are  stochastic  differential  equations 
with  time  random  coefficients.  In  the  latter  case  the  equations  of  motion 
are  differential  equations  with  random  constant  coefficients  [11]  or  with 
coefficients  that  vary  randomly  with  the  spatial  coordinates  (random 
fields).  Systems  with  parameter  uncertainties  are  referred  to  in  the 
literature  as  "disordered  systems."  The  methods  of  treating  dynamic 
systems  under  parametric  random  excitations  are  different  from  those  used 
in  solving  differential  equations  with  random  constant  coefficients  or  ran¬ 
dom  fields. 

In  problems  involving  random  constant  coefficients  the  engineer  is  con¬ 
cerned  with  three  main  problems:  random  eigenvalues,  random  response,  and 
optimum  design.  With  reference  to  linear  disordered  systems  these  three 
problems  will  be  reviewed  in  the  next  two  sections. 

IV . 2  Random  Eigenvalues 

The  eigenvalue  of  simple  single  degree-of-f reedom  systems  is  given  by 
the  square  root  of  the  stiffness  to  mass  ratio.  This  value  is  assumed 
constant  for  identical  systems.  However,  the  actual  eigenvalue  of  each 
system  deviates  from  the  original  calculated  one  because  in  reality  the 
physical  properties  of  the  elements  can  neither  be  measured  exactly  nor 


manufactured  exactly.  Thus,  the  eigenvalues  are  random  variables  whose 
statistical  properties  are  determined  by  the  random  coefficients  of  the 
inertia  and  stiffness  terms  of  the  equations  of  motion.  In  this  case  one 
may  be  interested  in  determining  the  probability  that  one  or  more  eigen¬ 
values  lie  in  a  given  range  or  less  than  a  certain  value.  Alternatively, 
one  may  need  to  know  the  probability  that  the  smallest  eigenvalue  is  spe¬ 
cified  in  a  given  range  [12].  Boyce  and  Goodwin  [13]  classified  parameter 
uncertainties  into  two  classes.  These  are: 

(i)  uncertainties  in  the  geometry  and  the  material  properties,  and 
(ii)  uncertainties  in  the  support  mechanism  of  the  system. 

These  uncertainties  appear  in  the  equations  of  motion  or  in  the  boundary 
conditions,  respectively.  Under  these  types  of  uncertainties  the  eigen 
value  has  been  determined  for  a  limited  class  of  dynamic  systems.  Boyce 
[12]  addressed  a  number  of  techniques  to  determine  the  statistics  of  the 
eigenvalues  of  systems  described  by  partial  differential  equations  and 
boundary  conditions  involving  uncertainty  in  their  parameters.  Two  mathe¬ 
matical  approaches  known  as  "honest”  and  "dishonest"  have  been  adopted  in 
the  literature  [14]  .  In  the  honest  approach  the  eigenvalues  are  first 
expressed  in  terms  of  the  system  parameters.  The  statistical  charac¬ 
teristics  of  this  solution  arc  then  determined  in  terras  of  the  statistical 
characteristics  of  the  random  parameters.  However,  this  approach  involves 
dif f  iculties  since  it  is  not  possible  to  express  the  eigenvalue  exactly 
except  for  very  few  simple  cases.  Four  honest  methods  are  outlined  by 
Boyce  [12].  These  are  the  variational  principles,  perturbation  methods, 
the  kernel  trace  estimates,  and  asymptotic  estimates.  In  the  dishonest 
approach  the  statistical  moments  of  the  eigenvalues  are  directly 
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determined  by  performing  averaging  analysis  to  the  system  partial  differen¬ 
tial  equations  and  its  associated  boundary  conditions. 

In  a  series  of  papers,  Purkert  and  Scheidt  [15-17]  established  a  number 
of  theorems  pertaining  to  functionals  of  weakly  correlated  processes. 

These  processes  are  encountered  in  the  eigenvalue  problems,  boundary  value 
problems  and  initial  value  problems.  Purkert  and  Scheidt  treated  the 
stochastic  eigenvalue  problem  for  ordinary  differential  equations  with 
deterministic  boundary  conditions.  The  coefficients  of  the  differential 
operator  were  assumed  to  be  independent  weakly  correlated  processes  of 
small  correlation  length.  As  the  correlation  length  vanishes  the  eigen¬ 
values  and  eigenvectors  were  found  to  possess  Gaussian  distributions.  In 
their  recent  monograph,  Scheidt  and  Purkert  [18]  treated  the  moments  of  the 
eigenvalues  and  mode  shapes  of  random  matrices  and  random  ordinary  dif¬ 
ferential  operators.  The  calculations  of  these  moments  were  based  on  per¬ 
turbation  expansions,  and  so  require  the  random  terms  to  be  appropriately 
small . 

Soong  and  Bogdanoff  [19]  examined  the  statistical  properties  of  the 
natural  frequencies  of  a  linear  n-degree-of-f reedora  system  whose  properties 
are  known  in  a  stochastic  sense.  They  used  a  method  based  on  the  transfer 
matrix,  developed  originally  by  Kerner  [20],  together  with  a  perturbation 
type  expansion.  For  a  linear  10  degrees-of-f reedom  system  with  random 
parameters  having  normal  distribution  with  small  standard  deviation  it  was 
found  that  the  top  few  natural  frequencies  have  values  which  are  very  sen¬ 
sitive  to  the  parameter  variations,  whereas  the  lowest  few  are  insensitive 
to  these  variations.  They  derived  explicit  expressions  for  the  natural 


frequencies  in  terms  of  the  parameter  variations.  These  expressions  can  be 
used  in  estimating  the  changes  in  frequencies  produced  by  deterministic 
parameter  changes.  Bliven  and  Soong  [21]  determined  the  statistics  of  the 
natural  frequencies  of  a  simply  supported  elastic  beam  with  random  imper¬ 
fections  in  the  beam  stiffness.  The  beam  was  modeled  as  a  lumped-parameter 
model  and  the  same  technique  of  Soong  and  Bogdanoff  was  employed.  Bliven 
and  Soong  found  that  when  the  stiffness  fluctuation  has  zero  correlation 
distance  the  natural  frequency  standard  deviation  vanishes.  The  standard 
deviation  reaches  the  value  of  0.5  when  stiffness  variation  is  perfectly 
correlated.  In  addition,  the  standard  deviation  of  the  beam  natural  fre¬ 
quency  was  found  to  be  insensitive  to  the  number  of  segments  in  the  lumped 
parameter  model. 

Collins  and  Thomson  [22]  treated  the  problem  of  eigenvalue  and  eigen¬ 
vector  statistics  of  a  simple  chain  of  equal  springs  and  masses  with 
uncorrelated  random  masses  or  with  random  uncorrelated  stiffnesses.  They 
showed  that  the  standard  deviation  of  the  frequency  is  governed  linearly 
with  the  standard  deviations  of  the  masses  and  stiffnesses.  This  result 
was  obtained  earlier  by  Soong  and  Bogdanoff  [19].  However,  these  linear 
relationships  disappear  when  a  correlation  exists  between  the  masses  and 
st i f  fnessps . 

Vaicai^i^  [23j  employed  a  two-variable  perturbation  expansion  procedure 
to  determine  the  eigenvalues  and  normal  modes  of  beams  with  random  and/or 
non-uniform  characteristics  which  do  not  deviate  considerably  from  the 
beam  mean  properties.  He  used  a  Monte  Carlo  simulation  to  determine  the 
statistical  averages  of  beam  eigenvalues  and  mode  shapes.  It  was  found  that 
the  eigenvalues  and  mode  shapes  deviate  significantly  from  those  of  a 
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uniform  beam.  The  difference  was  mainly  attributed  to  the  fact  that  gradual 
change  in  the  beam  stiffness  was  permitted.  In  this  case  the  beam  is 
"soft"  at  one  end  and  "hard"  at  the  other  end. 


IV. 3  Random  Response 

In  an  attempt  to  examine  certain  aspects  of  the  dynamic  behavior  of 
statistically  defined  systems,  Bogdanoff  and  Cbeanea  [24]  treated  linear 
single  degree-of-f reedom  systems  with  independent  discrete  distributions  in 
the  mass,  damping,  and  stiffness  coefficients.  Small  dispersion  in  the 
system  parameters  resulted  in  considerable  dispersion  in  the  system  fre¬ 
quency  response.  Their  analysis  was  based  on  a  partial  differential 
equation  for  the  response  joint  density  function.  This  equation  is  known 
as  the  Liouvi lie  equation  [11]  and  is  identical  to  the  Fokker-Planck 
equation  with  zero  diffusion  coefficient.  The  impulse  response  of  a  single 
dogree-of-f reedom  system  with  random  parameters  was  determined  by  Chen  and 
Soroka  [25]  by  using  a  perturbation  approach.  They  found  that  both  the 
mean  and  standard  deviation  of  the  response  were  non-stat ionary  and  the 
standard  deviation  was  90  degrees  out  of  phase  from  the  mean.  They 
concluded  that  for  systems  with  a  very  high  natural  frequency,  the  uncer¬ 
tainty  in  the  natural  frequency  has  a  very  negligible  effect  on  the 
response  statistics.  However,  the  effect  is  significant  if  the  natural 
frequency  is  low.  As  the  damping  factor  decreases,  the  dispersion  from  the 


mean  becomes  substantial.  In  another  study,  Chen  and  Soroka  [26]  con¬ 
sidered  the  response  of  multi-degree-of-f reedom  systems.  Their  study  indi¬ 
cated  that  the  response  statistics  of  disordered  systems  are  higher  than 


those  of  purely  deterministic  systems. 


V  V  ✓ 
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The  instantaneous  transient  response  statistics  of  an  undamped  linear 
multi-degree-of-f reedom  system,  subjected  to  arbitrary  but  deterministic 
forcing  functions,  with  stiffness  uncertainty  was  investigated  by 
Prasthofer  and  Beadle  [10].  For  the  case  of  an  impulsive  excitation  to  a 
single  degree-of-f reedom  system,  they  found  that  the  growth  of  the  response 
uncertainty  is  exponential.  As  the  standard  deviation  of  the  stiffness 
increases  the  response  mean  square  increases  rapidly  with  time.  For  a 
mult i-degree-of-f reedom  system  the  response  decay  ra.e  decreases  as  the 
correlation  coefficient  between  the  stiffness  elements  increases. 

The  influence  of  damping  uncertainty  on  the  frequency  response  of  a 
linear  multi-degree-of-f reedom  system  was  examined  by  Caravani  and  Thomson 
[27].  They  determined  the  mean  and  standard  deviation  of  the  response  by 
using  a  linearization  technique  and  a  Monte  Carlo  simulation.  They  pointed 
out  that  an  accurate  estimate  of  the  damping  coefficients  for  lightly 
damped  systems,  in  the  neighborhood  of  a  natural  frequency,  is  very  impor¬ 
tant  in  determining  the  mean  and  standard  deviation  of  the  system  response. 


IV. 4  Design  Optimization 

During  the  design  stage  of  structural  systems  the  fluctuations  of  their 
dynamic  characteristics  such  as  response  or  eigenvalues  should  be  defined. 
The  main  p’.uht  >ra  is  how  to  restrict  the  fluctuations  of  the  sysietn  para¬ 
meters.  For  example,  in  systems  in  which  the  values  of  displacement  are 
significant,  or  in  structures  for  which  the  safety  factors  for  fatigue 
strength  are  determined  in  terms  of  probability  functions,  the  problem  is 
to  set  up  an  optimum  standard  of  manufacturing  their  components.  Here  the 


permissible  fluctuation  in  the  characteristics  becomes  a  restrictive 


condition.  Under  these  circumstances,  the  designer  is  encountered  with  a 
problem  of  optimization  to  specify  the  maximum  permissible  fluctuations  of 
parameters . 

Tanaka  and  Onishi  [28]  developed  a  method  of  regulating  the  deviations 
of  random  parameters  and  derived  a  restrictive  conditional  formula  from  the 
permissible  fluctuation  of  displacements  or  natural  frequencies.  The 
method  is  based  on  the  linear  deviation  analysis  with  partial  differential 
analysis  together  with  sequential  linear  programming  (SLP)  for  a  number  of 
restrictive  conditions.  Tanaka,  et  al.  [29]  treated  the  optimization 
problem  of  the  allowable  variance  of  random  parameters  by  using  a  pertur¬ 
bation  method  and  Monte  Carlo  simulation.  They  computed  the  deviation  of 
the  steady  state  response  of  structural  systems  with  random  parameters  with 
the  purpose  of  regulating  the  deviation  of  the  random  parameters  when  the 
deviation  of  the  response  is  specified.  Rao  [30]  applied  the  multi¬ 
objective  optimization  techniques  to  the  design  of  simple  structures 
involving  uncertain  parameters  and  stochastic  processes.  The  necessity  of 
optimizing  the  structural  systems  involving  dynamic  restrictions,  random 
parameters,  stochastic  processes,  and  multiple  objectives  has  been  outlined 
by  Rao  [ 31 ] . 

In  this  chapter,  the  effects  of  randomly  t ime-varving  dampi  >t?  an; 
stiffness  (represented  by  stochastic  processes)  on  the  system  response  will 
be  investigated  by  using  the  Fokker-Planck  equation  approach.  Gaussian  and 
non-Gaussian  closure  schemes  will  be  used  to  obtain  the  mean  square 
response  of  the  system  in  the  normal  and  generalized  coordinates.  The  pre¬ 
sent  analysis  will  not  provide  any  information  about  the  stochastic  stabi¬ 
lity  of  the  system.  The  investigation  of  stochastic  stability  of 
aeroelastic  structures  will  be  examined  in  another  report. 
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IV. 5  Theoretical  Analysis 

Considering  again  the  aeroelastic  system  shown  in  Fig.  1.1,  and 
allowing  random  variations  in  the  damping  and  stiffness  coefficients, 
equations  (III. 3)  become 

Y^  +25^1  +  £c  (t)]y|  +  [l  +  gk  (t)]y1  =  -  ?  (x)a1 

-  f  [.j*!  +  a3V2]  +  +  a5Y,Y;  +  +  a^Y^  +  a/2 


+  VlY2  +  a10Y2 


(IV. 1) 


Y2  +  2c2r[l  +  Cc  (t)]Y2  +  r  [l  +  ^^^(t  )]  Y2  =  -  X  (t  )bt 


-  5  L  Vi  +  b3Y2J  +  VlYl  +  b5YlY2  +  b6Y2Yl  +  b7Y2Y2 

+  +  VlY2  +  b10Y22 

where  C  (x )  and  (t )  are  assumed  independent  Gaussian  random  processes 
C  .  K  , 

1  1 

in  damping  and  stiffness  in  the  ith  mode,  respectively,  (i  =  1,  2). 
Equations  (IV. 1)  can  be  approximated  through  successive  elimination  of 

from  non-linear  terms  and  will  be  transformed  into  the  Stratonovich  type 
equation  [1] 


X.  =  f.(X,x)  +  \  G.  (X,tX,(t  ) 
J  “ J 


(IV.  2) 


through  the  coordinate  transformation 


(IV. 3) 


equations  (IV. 1)  can  be  written  in  the  for  a 


X1  "  X3 


X2  "  X4 


(3  =  -  -  2^LX3  -  a4X3  -  (afi  +  r^X^  -  r2ayX2  -  %  ^X^ 

2c2ra5XlX4  -  2^^^  -  2?2rayX2X4  +  agX3  +  agX3X4  +  a1QX4 

-  Xl«k1(T)  "  *ltX3  +  a4XlX3  + 

~  r2[  a5X1X2  +  a7X2^k2(T)  '  2^2rta5XlX4  +  a7X2X4^c2(° 

-[a3  +  a2x3  +  a3x2  +  a4x2  +  a5xxx2  +  A6X2]5(t) 

(4  =  -  r2X2  -  2^rX4  -  b^X*  -  (  bg  +  r2b5)  X^  -  r2byX2  ( 

-  2clW3  -  %2tb5X1X4  -  -  2c2rb7X2X4  *  bjX2 

+  b9X3X4  b  »10X4  -  +  b6X2X3^c,(T) 


rz[x2  +  b5XtX2  +  b7X2]c^(T) 

2^2r[x4  +  b5XLX4  +  b^xj*  (t  } 

[Bt  +  B2X3  +  B3X2  +  B4X2  +  B5X1X2  +  Bf,X2jf,(i) 


Alternatively,  equations  (IV. 4)  can  be  written  in  the  Ito  form  [1] 


‘i  '  t  II  okJ 


9Gi  i 

II  Gi  •  dx  +  I  G.  . ( X ,x  ) d B .  (x  ) 

,  .  k  j  3  X,  1  h  1 1  ~  i 

kj  J  k  J 


n_\  <L7.  .--1 . ■ _ _ ^ 


where  Che  random  processes  ^  ^ (t )  have  been  replaced  by  white  noise  processes 

W^(x).  In  (IV. 5)  the  formal  derivative  of  the  Brownian  motion 
dBi(x) 

W.(t)  =  — ^ -  is  used.  dB^(x)  are  independent  Brownian  motion  processes 

with  properties 


E[dB.(x)]  =  0 

E[dB2(T>]  =  2D.  dr  =  o2dx 
L  1  J  l  i 


( IV .  6 ) 


Introducing  the  Ito  correction  term  represented  by  the  double  summation 
expression  in  equation  (IV. 5),  equatons  (IV. 4)  become 
dXt  =  X3dx 


dX2  =  X^dx 


dX3  =  {  -  xi  -  2?1X3  -  a4XL  *  (afi  +  r  a5)  X^  -  r2a?X^  -  2,  ^X^ 

-  2?2ra5X1X4  -  ^^X^  -  2c2ra?X2X4  +  agX2  +  agX3X4  +  a1QX2 

+  ^^c2^2r  [  a5XiX4  +  a7X2^4  a5B5XlX4  +  ^  a5^7  +  a7B5^XlX?X4 
+  a^b^X2X4]}dx  -  X^dB^  (x  )  -  r  [  a^X^X2  +  ayX2X4]  dB^  (x  ) 

-  2Mx3  +  a4X!X3  +  a6X2X3^dBCl(T) 

"  2rL  a5XlX4  H'  ^7X2X4]dBC2(x) 

A3  +  A2Xl  +  A3X2  +  A4X^  +  A5XlX9  +  A6X2]dB(x) 


(IV. 7) 
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dX,  =1-  r2X„  -  2z,  ~rX,  -  b,  X?  -  f  b,  +  r2bc)x,X„  -  r2b,X2 
4  1  2  24  41  ^6  5;  l  2  72 

-  2^1b4X1X3  -  2?2rb5X1X4  -  2^1b6X2X3-  2c2rb?X2X4 
+  bsx5  +  b,x3x4  +  bI0x4  +  3[  b^X^Xj  +  b6x2x3 

+  a4b4XlX3  +  (a4b6  +  V6)X1X2X3  +  a6b6X2X3] 

+  4D  r,  2r2f  X.  +  2bcX.X.  +  2b,X0X.  +  b2X2X.  +  2b.b,X,X„X. 
c2  2  L  4  514  724  514  57124 


+  ^Xjxjldr  -  2^[b4X1X3  +  b6X2X3]dB  (t) 

-  r2[X2  +  b5X1X2  +  b7X2]dBk^(T) 


-  2?2r[x4  +  b5XxX4  +  b7X2X4]dBc^(T) 


-  [Bx  +  B2X1  +  B3X2  +  B4X2  +  B5X3X2  +  B6X2]dB(x) 


(IV. 7) 
cont 'd . 


The  general  differential  equation  of  the  response  joint  moments  is 
obtained  by  using  the  Fokker-Planck  equation  approach  as  outlined  in  Chapter 

III.  This  procedure  results  in  the  moment  equation: 

1 

mi,j,k,!t  =  imi-l,  j,k+l,d  +  jmi,  j-l,k,*+l 

f  kt*  “l+l,J,k-l,*  +Cl^iDc  "  2j  mi,j,k,X  “ 


(a6  +  r  a5)rai+li  j+lik_1<£  ~  r  a7mi ,  j+2  ,k-l  ,£ 

+  241a4(c1D£;^-  1)tni+1)  j>k)£  +  ^2ra5^2rDc2"  2)  mi+l ,  j  ,k-l  ,£  +1 

+  2;Ta6(c  1DC^-  ljm^  j+L)k,x.  +  ^2ra7^2rDc2'  ^  rai  ,  j+1  ,k-l  ,£+1 

+  a8mi ,  j  ,k+l  ,1  +  a9rai ,  j  ,k,£+l  +  a10mi ,  j ,k-l  ,1 W 


1  i,j+l,k,£-l  2  c2  '  l,j,k,£  4  i+2,j,k,£-l 

-  (r  b5  +  b6)mt+1>  "  r  b7  rai ,  j+2  ,  k  ,£  - 1 

+  Clb4^1DCl"  2)mi+l,  j,k+l,£-l  +  2c2rb5^2rDc2"  ^“l+l.j.M 

+  Clb6^1DCl"  2K,  j+l,k+l,£-l  +  2c2rb7^2rDc1'  ^  mi ,  j+1 ,  k  ,Z 

b8mi,  j  ,k+2  ,£-l  b9mi,j,k+l,£  +  b10mi ,  j  ,k,£ +li 
+  k(k-n(DA2mij.ik_2^  +2DAlA2«i+lfj>k-2f)l 

+  2MlVi,j+l,k-2,£  + 

+  2D(AlA5+A2A3)ro.+1)  .+]^k_2^  +  A3+2A1 A6^  mi ,  j+2  ,k“2  ,£ 

+  4°Cl  A  -i.J.k.ti 

+  W{2DA1Blml(j>Ie+1>Jl.1  +  2d[  A^^bJ  m.+1  <  >£  ^ 

+  2D(A1B3+A3B1)m.ji+ljk_1>£_1 
+  2D(A1B4+A4BI+A2B2)m.+2>j>k_1>£_1 
+  2D(  At B5+A5Bl+A2B3+A3B2) m.+1 ^  j+1  ^  >£ ^ 
f  2D(  A3B3+AlB6+A6B1)  m.  ^  j+2  >k_y  >Jt  J 

+  Ut-1)l2DBK,j,k,I-2  +  2D»lB2»1+1,j,k,t.2 

+  2D(  B^+BjB-j)  m.+1  ^+1  tktl_2 

+  [\r2  +  +  40c2';2r2"i,.i,k.«> 


Again,  this  equation  reveals  that  the  response  moment  equations  constitute 
an  infinite  hierarchy  set.  This  infinite  hierarchy  may  be  closed  via 
Gaussian  or  non-Gaussian  schemes. 

For  the  Gaussian  closure  scheme,  fourteen  equations  of  first  order  and 
second  order  moments  will  be  generated  from  equation  (IV. 8).  These 
equations  are  coupled  with  third  order  moments,  and  if  the  third  order 
cumulant  is  set  to  zero,  the  fourteen  equations  can  be  closed.  The  response 
of  the  system  can  then  be  determined  by  integrating  the  fourteen  moment 
equations  using  the  IMSL  DVERK  routine. 

For  the  non-Gaussian  solution,  which  is  more  accurate,  69  moment 
equations  for  the  first  four  orders  are  generated.  These  equations  are 
found  to  be  coupled  through  fifth  order  moment  terms  which  can  be  closed  by 
setting  the  corresponding  cumulant  to  zero.  The  closed  69  moment  equations 
are  integrated  numerically  by  again  using  the  IMSL  DVERK  routine. 


IV. 6  Response  of  the  System  with  Damping  Uncertainty 


The  time  history  response  according  to  the  Gaussian  closure  is  plotted 

in  Fig.  IV. 1  for  internal  tuning  resonance  r  =  co 2 /ou ^  =  0.5,  damping  ratios 

~  ^2  ~  mass  ratio  n^/m^  =  0.2,  beam  length  ratio  1 2^  \  =  0.6,  and 

spectral  density  of  randomly  varying  damping  D  =  D  =  0.1  D  where  2D  = 

C1  C2 

0.08  is  the  spectral  density  of  the  base  motion.  It  is  seen  that  after  a 
long  period  of  the  time  parameter  x  =  1000  the  Gaussian  time  responses 
fluctuate  between  two  limits,  indicating  that  the  system  does  not  achieve  a 
steady  state. 

For  the  initial  conditions  indicated  on  each  figure,  the  first  mode 
grows  very  fast  with  rapid  fluctuations  while  the  second  mode  shoots  over 
its  initial  value  with  a  general  decay  which  reaches  values  below  the  ini¬ 
tial  value.  The  two  modes  exchange  energy  over  the  transient  period  t  = 
150.  During  the  steady  state  period  each  mode  fluctuates  between  two  enve¬ 
lopes.  The  period  of  oscillation  of  the  first  mode  envelope  is  twice  that 
of  the  second  mode.  Another  important  feature  is  that  the  level  of  the 
me;  n  squares  during  the  transient  period  is  higher  than  the  steady-state 
level . 

The  time  history  response  according  to  the  non-Gaussian  solution  is 
shown  in  Fig.  IV. 2.  This  figure  shows  fluctuations  during  the  i.^ansient 
response  period  r  =  160.  Contrary  to  the  Gaussian  solutions,  the  system 
response  achieves  a  stationary  state. 

The  effects  of  random  variation  in  danping  ratios  s  ^  and  c,  2  on  the 
Gaussian  and  non-Gaussian  autoparametr i u  region  in  normal  and  generalized 
coordinates  are  shown  in  Figs.  IV. 3  and  IV. 4  as  functions  of  the  internal 


6 

tuning  ratio  r  =  In  Gaussian  closure,  as  damping  is  varying  ran¬ 
domly,  the  autoparametr ic  interaction  region  becomes  more  narrow  and  the 
peak  values  of  the  two  mean  square  responses  stay  almost  the  same. 


Gaussian  closure  solution  for  system  with  randomly  varying  damping. 
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IV .7  System  Response  with  Stiffness  Uncertainty 

The  system  mean  square  responses  are  determined  in  this  section  when 

the  stiffnesses  of  the  system  are  varied  randomly  with  spectral  densities 

D,  =  D,  =  0.1  D,  (D  =  D  =0)  where  2D  =  0.08  is  the  spectral  den- 
\  K2  C1  C2 

sity  of  the  base  motion.  Figure  IV . 5  shows  the  time  history  responses 
according  to  the  Gaussian  closure  scheme  for  internal  resonance  r  =  0.5, 
damping  ratios  ^  \  =  ^  2  =  0.02t  mass  ratio  n^/m^  =  0.2,  beam  length  ratio 
i  2^1  =  0*6,  and  non-linear  coupling  parameter  e  =  0.02.  The  transient 
response  exhibits  the  same  features  outlined  in  section  IV. 6  and  the  steady 
state  response  fluctuates  between  two  quasi-stationary  envelopes. 

Figure  IV. 6  shows  the  time  history  responses  according  to  the 
non-Gaussian  closure  scheme  for  the  same  system  parameters  as  above. 

Again,  the  transient  response  fluctuates  with  energy  exchange  between  two 
modes  until  the  system  achieves  a  stationary  level. 

Mean  square  responses  obtained  by  Gaussian  and  non-Gaussian  solutions 
are  higher  in  the  transient  region  and  in  the  steady  state  than  the 
corresponding  values  when  the  system  possesses  damping  uncertainty  only. 
Figures  IV. 7  and  IV. 8  demonstrate  the  system  displacement  mean  squares  as 
functions  of  the  internal  tuning  parameter  r,  according  to  Gaussian  and 
non-Gens'-  an  solutions,  respectively.  It  is  seen  that  the  inclusion  of 
stiffness  uncertainty  has  a  remarkable  effect  on  the  mean  square  responses. 
The  mean  squares  of  the  response  displacements  are  relatively  higher  than 
those  without  stiffness  uncertainty.  In  addition,  although  the  region  of 
autoparametric  interaction  is  wider,  the  absorber  effect  is  less  pronounced. 
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tion  for  system  with  constant  damping  and  randomly  varying  stiffness 
damping  and  stiffness,  -------  randomly  varying  stiffness 


damping  and  randomly  varying  stiffness 


IV. 8  System  Response  with  Damping  and  Stiffness  Uncertainties 

In  this  section  the  combined  effect  of  the  system  damping  and  stiffness 

uncertainties  on  the  system  random  response  will  be  examined.  Typical 

spectral  densities  for  the  uncertainties  of  these  parameters  are  adopted, 

namely  D  =  D  =  D,  =  D,  =  0.1  D,  where  2D  =  0.08  is  the  spectral 
C1  C2  *1  K2 
density  of  the  base  random  motion. 

The  time  history  responses  according  to  the  Gaussian  closure  solution 
for  internal  tuning  ratio  r  =  0.5  are  shown  in  Fig.  IV. 9.  Corresponding 
plots  according  to  the  non-Gaussian  solution  are  shown  in  Fig.  IV. 10.  The 
inclusion  of  the  uncertainties  of  the  system  parameters  bring  the  system 
response  into  a  wider  fluctuating  limit  cycle  for  the  Gaussian  closure 
solution  or  into  a  wider  stationary  limit  cycle  for  the  non-Gaussian  clo¬ 
sure  solution.  Furthermore,  the  transient  response  period  with  parameter 
uncertainties  is  relatively  longer  than  the  transient  response  period  with 
constant  system  parameters.  The  mean  squares  of  the  response  displacements 
during  the  steady  state  regime  are  plotted  in  Figs.  IV. 11  and  IV. 12  by 
dotted  curves;  the  solid  curves  show  the  response  of  the  system  with 
constant  parameters.  The  higher  level  of  the  response  mean  squares  is  due 
mainly  to  stiffness  uncertainty  rather  than  damping  as  inferred  from  the 
previous  two  sections. 
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IV. 9  Conclusions 

The  influence  of  uncertainties  in  system  parameters  such  as  damping  and 
stiffness  on  the  autoparametric  response  of  an  aeroelastic  structure  has 
been  examined  numerically.  The  damping  influence  upon  the  response  charac¬ 
teristics  is  found  to  be  very  small  while  the  stiffness  uncertainty  has  a 
remarkable  effect  on  both  the  level  of  the  mean  square  response  and  the 
region  of  autoparametric  interaction.  The  transient  responses  in  the  time 
domain  are  relatively  higher  than  the  steady  state  responses.  This  feature 
is  very  important  in  determining  the  actual  stresses  under  random  vibra¬ 
tion.  With  stiffness  uncertainty,  the  system  response  uncertainty  becomes 
higher  than  when  the  system  possesses  constant  stiffness. 


CHAPTER  V 


CONCLUSIONS 

The  linear  and  non-linear  random  modal  interactions  of  a  two  degree-of- 
freedom  aeroelastic  structure  has  been  examined.  The  structure  was  sub¬ 
jected  to  a  random  wide  band  support  motion.  The  equations  of  motion  w _re 
derived  by  Lagrangian  formulation  and  transformed  into  normal  coordinates 
to  eliminate  the  linear  dynamic  coupling.  The  random  excitation  appeared 
as  parametric  and  non-homogeneous  terms  in  the  equations  of  motion. 

The  linear  modal  interaction  was  analyzed  via  the  Fokker-Planck 
equation  approach.  The  statistical  moments  of  the  system  response  were 
solved  for  various  system  parameters.  The  general  trend  of  the  linear 
solution  showed  that  there  is  a  suppression  of  one  mode  due  to  the  presence 
of  feedback  forces  from  the  other  mode.  Such  interaction  took  place  when 
the  frequencies  of  the  uncoupled  beams  are  close  to  each  other.  Further 
analysis  has  shown  that  the  parametric  terms  have  a  very  small  effect  on 
the  level  of  the  response  mean  squares.  However,  these  parametric  excita¬ 
tions  affect  the  stochastic  stability  of  the  system  equilibrium  con- 
f iguration. 

When  the  structure  is  tuned  to  the  appropriate  conditions  of  internal 
resonance,  the  non-linear  coupling  of  normal  modes  becomes  significant  in 
predicting  new  response  characteristics.  The  non-linear  analysis  is  not  a 
simple  task  and  involves  the  problem  of  infinite  hierarchy  of  differential 
equations  of  response  moments.  Two  schemes  have  been  used  to  close  the 
response  moment  equations.  The  first,  known  as  Gaussian  closure,  is  based 
on  the  assumption  that  the  response  does  not  depart  significantly  from 
Gaussian  distribution,  while  the  second  takes  into  consideration  the 


non-normality  of  the  response  and  is  known  as  the  non-Gaussian  closure 
method.  The  first  scheme  led  to  fourteen  differential  equations  in  the 
first  and  second  order  moments,  while  the  second  required  sixty-nine  dif¬ 
ferential  equations  in  the  first  four  orders  of  response  moments.  The  two 
sets  of  differential  equations  were  solved  by  numerical  integration.  Both 

solutions  exhibited  an  energy  exchange  between  the  two  modes  in  the  neigh- 

w2 

borhood  of  the  internal  resonance  condition  r  =  —  =0.5.  The  Gaussian 

W1 

closure  solution  gave  a  quasi-stationary  response  in  the  form  of  fluc¬ 
tuations  between  two  limits.  However,  the  non-Gaussian  solution  resulted 
in  a  strict  stationary  response.  The  accuracy  of  the  non-Gaussian  solution 
required  almost  twenty  times  more  computer  time  than  the  Gaussian  solution. 

In  Chapter  IV  the  damping  and  stiffness  coefficients  of  the  system  were 
subjected  to  time  random  variations.  It  was  found  that  the  damping 
variation  had  very  little  influence  on  the  response  characteristics,  while  the 
stiffness  variation  showed  a  pronounced  effect  on  the  response  mean  squares 
for  both  solutions. 

The  results  of  this  investigation  represent  part  of  an  ongoing  research 
project  supported  by  the  Air  Force  Office  of  Scientific  Research. 
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